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ABSTRACT 

We use the microlensing variability observed for eleven gravitationally lensed quasars 
to show that the accretion disk size at a rest-frame wavelength of 2500A is related to 
the black hole mass by log(iZ250o/cm) = (15.78±0.12)-h(0.80±0.17) log(MBH/lO9M0). 
This scaling is consistent with the expectation from thin disk theory (i? oc M^'^), but 
when interpreted in terms of the standard thin disk model (T oc R^^^^), it implies 
that black holes radiate with very low efficiency, log(r/) = —1.77 it 0.29 + \og{L/LE) 
where r/ = L/{M(?). Only by making the maximum reasonable shifts in the average 
inclination, Eddington factors and black hole masses can we raise the efficiency estimate 
to be marginally consistent with typical efficiency estimates (7/ ~ 10%). With one 
exception, these sizes are larger by a factor of ~ 4 than the size needed to produce 
the observed 0.8 //m quasar flux by thermal radiation from a thin disk with the same 
T oc temperature profile. While scattering a significant fraction of the disk 

emission on large scales or including a large fraction of contaminating line emission can 
reduce the size discrepancy, resolving it also appears to require that accretion disks have 
fiatter temperature/surface brightness profiles. 

Subject headings: accretion, accretion disks — dark matter — gravitational lensing: 
micro — gravitational lensing: strong — quasars: general 
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Introduction 



Despite nearly 40 years of work on accretio n disk physics, the simple IShak ura & Sunyaey 



(119731) thin disk model, its relativ istic cousins (e.g. iPage Sz Thorne 



1974 



1997 



Hubenv fc Hubenv .19971: 



Hubenv et al.ll2001 : Li et al.ll2005 ) and more sophistic ated implementations (e.g. lNaravan. Barret Sz McClintock 



De Villiers. Hawlev &: Krolik 



servational reservations (see lFrancis et al 



2003 



Blaes 



1991 



20071) remain the standard model despite some ob 



Koratkar fc Blaeslll999l : IColhn et al.ll2002l ). Quasar 



accretion disks cannot be spatially resolved with o rdinary telescopes, so we have been forced to 



test accretion physics through ti me variability (e.g. IVanden Berk et al.ll2004l: Sergeev 



2002^ : lBonning et al ] l2007l l. One notable success is the use of reverberation mapping (e.g 



et al 



2005 



Cackett. Horne &: Winkled 120071 ) and spectral modeling (e.g. ISun fc MalkanI Il989l : ICollin et al 



Peterson et al 



20041 ) of quasar broad line emission to calibrate the relation between emission line widths and black 



hole masses. The line emission, though, comes from scales much larger than the accretion disk, and 
attempts to use similar methods on the continuum emission have had limited s uccess, largely be- 
cause quasars show little optical variability on the disk light-crossing timescale (jCollin et al.ll2002l : 
Sergeev et aPbood ). 



Gravitational telescopes do, however, provide the necessary resolution to study the structure 
of the quasar continuum source. Each gravitationally lensed quasar image is observed through a 
magnifying screen created by the stars in the lens galaxy. Sources that are smaller than the Einstein 
radius of the stars, typically ~ 10^^ c m, show time varia ble fluxes whose amplitude is determined 
by the source size (see the review by IWambsganssI 120061 ) . Smaller sources have larger variability 
amplitudes than larger sources. In this paper, we exploit the optical microlensing variability ob- 
served in eleven gravitationally lensed quasar systems to measure the size of their accretion disks, 
and we find that disk sizes are strongly correlated with the masses of their central black holes. 

In ^we describe the monitoring data, the lens models we use based on Hubble Space Telescope 
(HST) images of each system and our microlensing analysis method. In ^we describe our accretion 
disk model and our results for the relationship between disk size and black hole mass. While we 
analyze and discuss the results in terms of a simplified thin disk model, they can be compared to 
any other model by com paring our measu r ement of the half light radius to that expected from the 
model of choice because Mortonson et al. (j2005l ) demonstrate that the half-light radius measured 
from microlensing is essentially independent of the assumed disk surface brightness profile. The 



^Based on observations obtained with the Small and Moderate Aperture Research Telescope System (SMARTS) 
1.3m, which is operated by the SMARTS Consortium, the Apache Point Observatory 3.5-meter telescope, which 
is owned and operated by the Astrophysical Research Consortium, the WIYN Observatory which is owned and 
operated by the University of Wisconsin, Indiana University, Yale University and the National Optical Astronomy 
Observatories (NOAO), the 6.5m Magellan Baade telescope, which is a collaboration between the observatories of 
the Carnegie Institution of Washington (OCIW) , University of Arizona, Harvard University, University of Michigan, 
and Massachusetts Institute of Technology, and observations made with the NASA/ESA Hubble Space Telescope for 
program HST-GO-9744 of the Space Telescope Science Institute, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under NASA contract NAS 5-26555. 
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surface brightness p rofile is best probed by measuring the dependence of the micro l ensing amplitude 
on wavelength (see Unguita et al.1 boosl: IPoindexter et al.1 bood : Ugol et alj I2OO9I : Bate et al.1 bood : 
Flovd et alj l2009l : iMosquera et al.1 12009| ) . In ^ we discuss the results and their implications for 
thin accretion disk theory. All calculations in this paper assume a flat ACDM cosmology with 
h = 0.7, flM = 0.3 and = 0.7. 



2. Data and Analysis 

We collected new monitoring data in the systems HE 0435-1223, FBQ 0951+2635, HE 1104- 
1805, SDSS 1138+0314, SBS 1520+530 and Q 2237+030 i n the r-, R- and F-bands on the 
SMARTS 1.3m using the ANDICAM optical/infrared camera (jPePov et ahlbood n. the Wisconsin- 



Yale-Indiana (WIYN) observatory using the WIYN Tip-Tilt Module (WTTMjfl, the 2.4m Hilt- 
ner telescope at the MDM Observatory using the MDM Eight-lfl, Echelle and RETRQCAIvi 



([Morgan et al.ll2005l ) imagers and the 6.5m Magellan Baade telescope using IMACS ( Bigelow et al 



I999I') . We supplemented our monitoring data with publis h ed quasar light curve s from lParaficz et al 



2000al ]bh. lGavnuUina et al.1 toom . and 



(I2OO6 ). Schechter et al.1 (ll997l'l.lWvrzvkowski et al.N2003 ) lOfek fc Maod ((20031), OGLE (jWozniak et al 



Poindexter et al 



(I200J. We measured the flux of each im- 



age by comparison to the flux from refer ence stars in the field of each lens. Our analysis of the 
monitoring data is described in detail by iKochanek et al.l ([200a) • In systems with published time 
delays, we offset the light curves by the delays to eliminate the intrinsic source variability. In this 
paper we also r nake use of previousl y reported optical ac cretion disk sizes for the lensed quasars 
QJ 0158-4325 dMorgan et al.l boOSal'l. SD SS 0924+0219 dMorgan et al.1 [200^1. SDSS 1 004+4112 



(IFohlmeister et al.1 \2008i ). PG 1115+080 (|Morgan et al.1 l2008al ^ and RXJ 1131-1231 (|Dai et al 



2OIOI). SDSS 0924+0219, SDSS 1138+0314 and Q 2237+030 do not have pubhshed time delays 
but are all quadruply lensed quasars with short (< 15 days) estimated delays. Since the typical 
timescale for microlensing is significantly longer than this (> 1 year), we ignored the delays in 
these systems and used their raw lightcurves in our microlensing analysis. For QJ 0158-4325 we 
developed an analysis method that allowed us the simultaneously estimate the tir ne delays and 
disk siz es including their mutual uncertainties, the details of which are described in lMorgan et al.l 
(l2008al l. 



All eleven lenses have been observed in the V- (F555W), /- (F814W) and ff-bands (F160W) 
using the WFPC2, ACS/WFC and NICMOS instruments on HST. We fit these images as com- 
binations of point sources for the quasars and (generally) de Vaucouleurs models for the lenses 



^http://www.astronomy.ohio-state.edu/ANDICAM/ 
■^http: //www. wiyn.org/wttm/WTTM_manual. html 
^littp://www. astro. columbia.edu/ arlin/MDM8K/ 
■■^littp://www.astronomy.ohio-state.edu/MDM/RETROCAM 
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as described in iLehar et alj (j2000l ) . These provided the astrometry used for lens models and de- 
fined a constant mass-to-light (M/L) ratio model for the mass d istribution iii the lens models. 
We modeled each system using the GRA VLENS software package (jKeetonl l200ll ) . For all systems 
except the cluster lens SDSS 1004+4112, we generate a series of ten mode ls starting from a con- 
stant M/L model and then adding an NFW ([Navarro. Frenk &: WhitdllQQq ) halo. The sequence is 
parametrized by /m/Lj the mass fraction represented by the visible lens galaxy relative to a con- 
stant M/L model. In general, we start with the constant M/L model, Jm/l = then reduce 
its mass in increments of A/^^/j;^ = 0.1 with the NFW halo's mass risin g to compensate. For th e 
cluster-lensed quasar SDSS 1004+4112, we use the fixed mass model from lFohlmeister et al.l (|2007l ). 
and we assume a set of 10 evenly spaced stellar mass fractions in the range 0.1 < Av*/k < 1.0. Thus, 
our results marginalize over any uncertainties in the dark matter halos of the lenses. For typical 
lens es, which we expect to be dark niatter dominated, microlensing favors low stel l ar mass fractions 



f e.g. Morgan et al.ll2008bl : iDai et al.ll2010l : IChartas et al 



200S 



Poolev et al. 



200S 



Mediavilla et al 



20091 ). We adopt, however, a very conservative philosophy and include a broad spectrum of models 
even though a stronger prior favoring dark-matter dominated models is relatively easy to justify. 
We do this so that our estimates of the source sizes depend as little as possible on other model 
assumptions. We have nonetheless tested the sensitivity of the size estimates on lens model dark 
matter fraction priors, and we have found very weak, if any, dependence. 

These lens models provide the convergence /t, shear 7 and stellar surface density needed 
to define the microlensing magnification patterns. We assume a lens galaxy stellar mass function 
dN{M)/dM (X M~i-3 with a dynamic range of a factor of 50 that ap proximates the Galactic disk 
mass function of Gould (2000, see also lPoindexter &: Kochanekll2010al ). We also know from previous 
theoretical studies that the choice of the mass function will have little effect on our conclusions 



given the other sources of uncertainty (e.g. iPaczynskil Il986l : IWyithe et al.ll200d ). For the typical 
lens we generated 4 magnification patterns for each image in each of the 10 lens models. We gave 
the magnification patterns an outer scale of 20{Re), where {Re) is the Einstein radius for the mean 
stellar mass (M). This outer dimension is large enough to fairly sample the magnification pattern, 
while the pixel scale of the 4096^ magnification patterns is small enough to resolve the accretion 
disk. We determined the properties o f the accret i on dis k by mod eling the observed ligh t curves using 
the Bayesian Monte Carlo method of lKochanekl (|2004l ) (also see lKochanek et al.ll2007l ). For a given 
disk model, we randomly generate light curves, fit them to the observations and then use Bayesian 
methods to compute probability distributions for the disk size averaged over the lens models, the 
likely velocities of th e observer, lens, s ource, stars and the mean microlens mass (M). We used the 
velocity model from iKochanekl (j2004l ). which used the projected CMB dipole (jKogut et al.lll993l ) 
for the observer, a stellar velocity dispersion set by the Einstein radius of the lens and peculiar 
velocity scales for the lens and source of 235/(1 + z) km s~^. We use a prior on the mean microlens 
mass of 0.1M(T^ < (M) < l.OMc^, but the disk size estimates are insensitive to this assumption (see 
Kochaned I2OO4I : iMorgan et aPboOSbl ^ . 



An obstacle to our analysis technique is how to handle differences between the observed flux 
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ratios and the flux ratios of the trial lightcurves. In systei ns where we are confident that the 
intrinsic flux ratios are known (e.g. PG 1115+080, for which IChiba et al.ll2005l . measured mid-IR 
Hux ratios), we set a prior to favor model lightcurves within 0.1 mag of the intrinsic flux ratios. In 
most cases, however, the intrinsic Hux ratios are not known because all observed Hux ratios can be 
affected by chromatic microlensing or dust extinction. In these cases we set a much looser prior 
on this magnification offset, typically ~ 0.5 mag. Loosening the prior on the magnification offset 
degrades our ability to to estimate the d ark matter conte nt of the lens galaxy, but it has little effect 
on accretion disk size measurements (see lDai et al.ll2010l . for a detailed demonstration of this effect 
in RXJ 1131-1231). 

We use black hole mass estimates for the quasars that are based on observed quasar emission 
line widths and the locally c alibrated viri a l relat ions f or black hole rnasses , for which we adopt the 
combined normalizations of Onken et al. ( 2004 ) and Greene &: Ho (j2007l ). For most systems we 
simply used the black hole mass estimates from IPeng et al.l ( 20061 ) based on the C IV (A1549A), 
Mg II (A2798A) and H/3 (A486lA) mass-linewidth relations. For SDSS 1138+0314, we mea- 
sured the width of the C IV (1 549A) line in optical spectra from the Sloan Digital Sky Survey 
dAdelman-McCarthv et ahlbood ) and for Q 2237+030 we used the C IV line width measurement 



from 



virial relation of 



Yee h De Robertid (Il99ll). We estimated the black hole masses for these systems using the 



Vestereaard fc PetersonI (120061). For SDSS 1004+4 112, we measured the wid th of 



the Mg I I fA2798A) emission line in spectra from llnada et al.l (j2003l ) and lRichards et al.l (j2004l ) and 
used the iMcLure &: Jarvia (|2002l ) Mg II virial relatio n to estimate its black hole mass. These mass 
estimates are reliable to approximately . 3 dex (see iMcLure &: JarvisI 120021 : iKollmeier et al.l l2006l : 
Vestergaard k. PetersonI l2006l : IPeng et al.ll2006l ). 



3. Results 



We model the surface brightness profile of t he accretion disk as a po wer law temperature 
profile, T oc matching the outer regions of a lShakura &: SunvaevI ( 19731 ) thin disk model. We 

neglect the central depression of the temperature due to the inner edge of the disk and corrections 
from general relativity to avoid extra parameters. The effect of this simplification on our size 
estimates is small compared to our measurement uncertainties provided the disk size we obtain is 
significantly larger than the radius of the inner disk edge. We will compare three disk size estimates 
in the context of this simple model. First, there is our size measurement from the microlensing, 
Rs- This microlensing size should be viewed as a measurement of the half-light radius, but we 
parametrize the results in terms of the simple thin disk model in order to facilitate comparisons to 
the thin disk model. Converte d to a half-light radiu s, i?i/2 = 2.44i?s, the measurements will be 
nearly model-independent (see iMortonson et al.ll2005l ). Second, there is the theoretically expected 
size as a function of black hole mass in the thin disk model (see Eqn. [2] below). Third, there 
is the thin-disk size which would yield the observed optical fiux assuming thermal radiation and 
a T oc R~^/^ temperature profile (see Eqn. [3] below). We emphasize that while the second and 
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third size estimates are based upon the same accretion disk model, they diff er in that th e y take 



Cohin et al 



different measured quantities, either Mbh or /-band ffux, as an argument 
had previous l y not iced that the theory size (Eqn. [2|) and the flux size (Eqn. [3]) may be discrepant. 
Poolev et al.l (|2007l ) argued that the microlensing sizes and the theory sizes (Eqn.[2|) are discrepant, 
but they used black hole masses based largely on estimates of the bolometric luminosity, which in 
some sense forced a reconciliation of the flux and theory sizes. 



We assume that the disk radiates as a black body, so the surface brightness at rest wavelength 



^rest is 



where the scale length 



_ 2hpC 

T3 

^rest 



exp 



R 



3/4 



(1) 



R 



45GXf^,tMBHM 



^rest 

167r6/ipc2 



1/3 



9.7 X 10^^ 



Vest 



/xm 



4/3 



/ Mbh 
V1O9M0 



2/3 



L 

TjLl 



1/3 



cm 



(2) 



is the radius at which the disk temperature matches the wavelength, kT\^^^_^ = hpc/\.est^ hp is 
the Planck constant, k is the Boltzmann constant, Mbh is the black hole mass, M is the mass 
accretion rate, L/Le is the luminosity in units of the Eddington luminosity, and r/ = L/{M(P') is 
the accretion efficiency. We can also compute the size under the same model assumptions based on 
the magnification-corrected I-band quasar fluxes measured in HST observations as 



Ri = 2.83 X 10 



15 



1 



COS I 



Dr 



Os \ f ^I,obs \^^^ -[^Q~0.2(J-19) ^-1 



cm 



(3) 



where Dos/th is the angular diameter distance to the quasar in units of the Hubble radius, / is 
the magnification-corrected magnitude and i is the disk inclination angle. 

Our results are shown in Figures [1] through [3] and summarized in Table [1] For the comparison 
with theory and the figures, we corrected the measured sizes to Xrest = 2500A assuming the X^^^ 
scaling of thin disk theory and a mean inclination (cosz) = 1/2. We chose 2500A because it 
was typical of the actual rest-frame wavelength (see Table 1), minimizing the sensitivity of our 
estimates to any uncertainty in the true wavelength scaling. Only the size of RXJl 131-1231 is 
strongly affected by changing the scaling of size with wavelength, because of its remarkably low 
source redshift ( z,^ = 0.658). Howev er, even if we vary the temperature profile of the disk from 

) to T oc R~^, corresponding to the range from Rx oc to 



le.g. 



Francis et al 



1991 



T oc , 

Rx oc A, the wavelength-corrected disk size changes only by ~ 25% and the fit parameters for the 
relationship between disk size and black hole mass change by less than 5%. 

Although we use a face-on disk model, we need to consider the role of inclination. Both the 
microlensing size and the flux size (Eqn. [3]) are set by the projected area of the disk, so the true 
microlensing disk scale should be 1 /V cos i of the face-on estimate. When we average over an 
ensemble of systems assuming a random distribution of inclination angles, we are averaging over 
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projected areas (i?^ 



true 



SO we correct our measurements to 



'fit/ - \-HrueCOSi) = (1/2 ± l/2y/3)Rl^^, 

{Rtrue) = V^Rfit- We discuss the consequences of other incHnation distributions in §|31 In Tabled] 
and Figs. [1] and O we use this average correction for both the microlensing and flux sizes. The gray 
band in Fig. [1] shows the expected extra variance of 1/2 V3 arising from the inchnation if we view 
our fits as matching the predicted and observed projected areas. 

There are two striking facts illustrated by the figures. First, we clearly see from Fig. [T]that 
the microlensing sizes are well correlated with the black hole mass. A power-law fit between -R2500 
and Mbh including the uncertainties in both quantities yields 



log 



^2500 

cm 



^15.78 ± 0.12) + (0.80 ± 0.17) log 



M, 



llO^M, 



BH 





2/3 N 



(4) 



which is consistent with the predicted slope from thin disk theory {R oc M^^). If we fix the slope 
with mass to 2/3 (see Fig. th e relation implies a typical Eddington factor of log{L/r]LE) = 
1.77 lb 0.29. iKollmeier et al.l (|2006l ) estimate that the typical quasar has L/Le ~ 1/3, which would 
indicate a radiative e fficiency of r? = L/^Mc?) ~ 0.006. This efficiency is very low compared to 



standard models (e.g. lGammidll999l ) and observational constraints on radiative efficiency derived 
from the local black hole mass and quasar luminosity functions (e.g. lYu &: Tremaind l2002l : ISoltan 
19821 ) . We discuss this point further in § [H 



The goodness of fit, = 9.15 for 9 degrees of freedom {x^ /N^of = 1.07) without including 
any effects from the spread in inclination, suggests that the errors on the size and Mbh esti- 
mates are appropriate. While the formal uncertainties in Mbh from the line width relations 
are only ^0.1 dex, th e systematic uncertainties are generally believed to be closer to 0.33 dex 



(|McLure k Jarvisll2002l : IKollmeier et al.ll2006l : IVestergaard k Peterson! 12009 : IPeng et al.1 12OO6I I . If 

we fit the relation assuming that this is an additional random error in each black hole mass, 
then we find log (i?250o/cm) = (15.9 ± 0.2) + (1.0 ± 0.3) log (Mb^/IO^Mq) with a goodness of fit 
/Ndof = 0.59. This implies either that 0.33 dex of random uncertainty for each individual black 
hole mass is too large or that we have over-estimated the uncertainties in the size measurements. 



The low efficiency estimate is closely related to the argument by iPoolev et al.l (|2007l ) that 
the sizes estimated from microle nsing are lar g er th an expected from thin disk theory, but the 
comparison is not exact because iPoolev et al.l (j2007l ) used black hole masses determined largely 
from estimates of the bolometric luminosity rather than the emission line width method. This 
causes some covariance between th eir theory and flux s izes so an exact comparison to our results is 
difficult. Nonetheless, compared to Pooley et al. (j2007l ) we find a smaller discrepancy because our 
full calculations of the microlensing sizes ten d to be somewhat s maller (by an average of 0.2 dex 
for our four common systems) than those of iPooley et al.l (120071 ). Despite these differences, our 
low estimate for the radiative efficiency is essentially the same as the problem pointed out by these 
authors. 



Second, there is a much more striking discrepancy between both the microlensing and theory 
size estimates (Eqn. [2]) with the flux sizes (Eqn. [3]). While most of the offset between the mi- 
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crolensing size measurements and the expectations from thin disk theory could be explained by the 
existing uncertainties, the offset from the estimate based on the quasar flux is more significant. The 
measured disk sizes are significantly larger than the flux sizes in all systems except QJ 0158-4325 
(see Fig. [3]), with an average offset of 0.6 ± 0.3 dex. Simply put, the quasars are not sufficiently 
luminous to have the sizes estimated from microlensing while radiating as black bodies with a 
T oc R^^/'^ temperature profile. 



Discussion 



Using microlensing we have made the first observational demonstration of the dependence of 
accretion disk sizes on black hole mass. While the slope of the relation is still relatively uncertain, 
it is consistent with the scaling of expected for quasars with similar Eddington ratios. The 

absolute scales correspond to relatively low radiative efficiencies, which we discuss in detail below. 



Like Pooley et al. (j2007l ) we find that the microlensing estimates of the disk size are somewhat 
larg er than would be exp ected from thin disk theory for typical Eddington ratios {L/Le ~ 1/3, 



e.g. 



KoUmeier et al.ll2006l ). We can examine this in terms of our estimate of the radiative efficiency 



log = (-2.25 ± 0.29) + log (^3-^^ + 2 log 



+ -log(2(cos i)), 



(5) 



where 3L/Le is the Eddington ratio normalized to 1/3, M b h, true I^Ibh, est is the average ratio of 
the true black hole masses to our estimates, and 2(cosz) is the dependence on changes in the mean 
inclination angle from that for a uniform distribution. Constraints de rived from inatching the total 
radia ted energy of quasars to the local black hole mass function (e.g. ISoltanlll982l : lYu &: Tremaine 
2OO2I ) argue for a typical radiative efficiency close to 10% or logry = — 1 and even higher efficiencies 
for the most luminous qu asars. It is bare l y possible to reconcile our estimates with these. First, 
in unification models (e.g. lAntonuccilll993l : lUrrv fc Padovanilll995l ) the inclination angles of optical 
quasars are preferentially face-on rather than uniformly distributed. Inde ed, the first microlensing 
measu rements of a disk inclination angle favor face-on viewing angles (jPoindexter &: Kochanek 
2010bl ). If, for example, we assumed quasars were uniformly distributed only over the range 1/2 < 



cosi < 1, then our estimate of the efficiency rises by 0.2 dex. The methods used in lKoUmeier et al 



(|200d ). or any other study of Eddington factors, will generally be more reliable for the distribution 
of the Eddington factors than for the absolute values. The absolute values are sensitive to the 
absolute calibration of the black hole masses and correctly estimating the bolometric luminosity. 
Shifting the typical quasar from L/Le = 1/3 to L/Le — 1 would reduce the discrepancy by 0.5 dex. 
Similarly, most studies of estimating black hole masses from emission line widths would accept that 
there are absolute calibration uncertainties of order 0.3 dex in MBH,true/MBH,est- Unfortunately, 
we gain only this factor rather than its square if we make MBH,true/MBH,est = 3, despite Eqn. \5\ 
because the estimates of the Eddington factor also have to be adjusted downwards if we raise the 
estimated black hole mass. If we move all three of these terms in the same direction {L/Le = 1 
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instead of 1/3, M bh, true /Mbh, est = 3 instead of 1 and (cosi) 
bring logr/ = —1.29 it 0.29 within Icr of 10% efficiency. 



3/4 instead of 1/2) we can we 



This may be stretching the allowed parameter shifts, but it does not address the still larger 
discrepancies between either the microlensing or the theory sizes with the flux sizes. To emphasize 
this point, if we estimate from the flux sizes we find logry = —0.43 ± 0.21 which has problems of 
the opposite sign. Are these discrepancies between these three size estimates due to a problem in 
the measurements, an oversimplification of the disk model or a fundamental problem in the thin 
disk model? We have tested our approach using Monte Carlo simulations of light curves and verified 
that we recover the input disk sizes. Our results are also only weakly sensitive to the assu med prior 



on th e microlens masses (see iKochanekl l2004l : [Morgan et al.l l2008bl . for a discussion) . iDai et al 



(|2010l ) in their detailed study of microlensing in RXJ 1131—1231 show that none of the details of 
our approach significantly affect the size estimates. This suggests that we must look to significant 
changes in the physical model to explain the differences. Fig. |4] illustrates the consequences of four 
possible modifications: scattered light, contaminating emission, effects of the inner disk edge, and 
changes in the temperature profile. In each case we rescale the ratio of the fiux and microlensing 
sizes while holding the optical fiux and the half light radius of the disk fixed to see if we can shift 
the ratio from the observed —0.6 ± 0.3 dex. Holding the half-light r adius fixed should mean that 
the new model will be consistent with the microlensing constraints (jMortonson et al.ll2005l ). We 
generalize the temperature profile of Eqn. [T]to 

1/4 

1 



r(K)«|f) 



/3 



(6) 



with locally thermal emission for R > Rm- Our fiducial models have /? = 3/4 and Rin = 0. Here 
Rx is just a scale length and does not correspond to the expression in Eqn. [2j Because Eqn. [6] is 
not based on a theoretical model, we have no means of relating Rx to physical parameters. We can 
attempt to reconcile the fiux and microlensing sizes with Eqn. [6l but we cannot use it to examine 
the efficiency problem. 

We have assumed so far that the observed optical flux is emission directly from the accretion 
disk, hence we will overestimate the microlensing size of the disk if an appreciable fraction of the 
emission originates on scales larger than the accretion disk. For example, if we flt the data modeling 
30% of the observed light as unmicrolensed emission from scales much larger tha n the a ccretion disk, 
then we flnd that the microlensing size estimates shrink by 20-50%. iDai et al.l (|2010l ) investigated 
this problem in detail for RXJ 1131—1231. As the fraction of contaminating light increases, the disk 
has to become more compact in order to produce the same amplitude of microlensing variability. 
This large scale emission can be either scattered emission from the disk or contaminating emission 
from some other source. The two effects differ because the flux from the disk is unchanged by 
scattering, while the flux from the disk is reduced by the amount of contaminating emission. 

The main sources of c ontamination will be emission from t he larger and minimally microlensed 
line emitting regions (e.g. ISluse et al.l 120071 : ISugai et al.l 120071 ) and the unmicrolensed emission of 
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the quasar host galaxy. Here we need not worry about the host galaxy as we are examining the 
rest frame ultraviolet emission of luminous quasars. The sources of line contamination include 
not only the obvious broad /narrow lines but also the broad Fe II and Balmer continuum emission 



(jMaoz et al.lll993l : IVestergaard Peterson 
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5I) thai 



c an represent 



^ ^ . 30% of the apparent con- 

tinuum flux at some wavelengths (|Netzer fcW illis IQsi: brandilllQsd l. Table □ notes the possible 



sources of broad line contamination for each system. While four of the systems have the Mg II line 
in the filter band pass, the Balmer continuum and Iron line complexes are probably the dominant 
source of line contamination. However, the line emission is generally reprocessed harder radiation, 
so as we increase the line contamination to reduce the microlensing size of the disk, we also reduce 
the optical flux coming from the disk and hence the flux size of the disk. Thus, producing the large 
scale emission by scattering the disk emission reduces the microlensing size but leaves the flux size 
unchanged, while producing it by line contamination reduces both size estimates. Fig. 0] shows that 
in our simple model, contamination and scattering cannot bring the two sizes into agreement unless 
most of the emission is not coming directly from the disk. Our discussion of scattering here assumes 
it is not accompanied by significant changes in the photon energy. If much of the observed optical 
radiation is due to Compton scattering of softer photons, then the net effect would depend on the 
physical scale of the scattering medium, particularly since the highest densities of hot electrons 
are likely to be on similar scales to that of the accretion disk. One consequence of reducing the 
directly observed emission from the disk is that the fractional variability of the disk emission rises 
in proportion, as the (re)emission on large (parsec) scales cannot produce the short time scales of 
the intrinsic quasar variability. This consideration probably rules out models in which a minority 
of the observed flux comes directly from the disk. 

The problem is also not a consequence of the simplifications in the disk model, namely the 
neglect of the inner edge and the different inner temperature profile of a relativistic disk (e.g. 
Page &: Thorndll974l ). Whether we use the microlensing or flux scales for the disk, the optical 
emission is mostly radiated fa r from the scale expected for the inner edges of disk (~ GOr^ rather 
than a few = GMbh/c^)- iDai et al.l (|20ld ) examined this problem for RXJ1131— 1231 using 
the full relativistic iHubenv et al.l ((20011) models and found few changes from our simple standard 
model. We illustrate this here by adding an inner edge to the disk with Rin = 0.1R\, which is 
larger than we would expect from the sizes and black hole masses of these sources. As shown in 
Fig. m this has little effect on the size ratio unless the disk temperature profile is very steep. 

The final possibility we consider is changing the temperature profile of the disk. As shown in 
Fig. m using a temperature profile flatter than /3 = 3/4 has the strongest effect on the size ratio 
of all the changes we consider. A temperature slope closer to /3 ~ 0.5 can bring the two scales 
into agreement. With the addition of a scattered or contaminating component, smaller changes are 
needed, and the inner disk edge becomes less important for the flatter profiles. Steeper temperature 
profiles, on the other hand, worsen the problem, although in this regime the inner disk edge also 
becomes important. 



Fig. Ulalso shows the existing microlensing limits on (3 from lAnguita et al.l (|2008l ). iBate et al 
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(12003), lEigenbrod et al.1 ^2008 ). iFlovd et alJ (120091) . and IPoindexter et al.l (|2008l ). based on the 
wavelength dependence of the microlensi ng. Mosauera et al. ( 20091 ) also report results consistent 
with /3 ~ 0.75. With the exception of iFlovd et al.l (j2009l ). these initial studies are consistent 
with both /3 = 3/4 and the shallower p rofiles that would help to resolve the differences in the 
size estimates. The discrepant result of iFlovd et al.l (j2009l ) r equiring a s teeper slop e is likely a 



conseq uence of their approach. Unlike the other studies, both iBate et al.l (120081 ) and iFloyd et al 



([20091) use only the color differences between images observed at single epochs rather than light 
curves. This means that at any wavelength they can only set an upper bound on the source size - 
it is the amplitude observed during caustic crossings that sets the lower bounds, and that can only 
be measured using time variability. However, given only a set of upper limits on the sizes that are 
smaller for shorter wavelengths, the estimate for /? will be det e rmine d by the priors used fo r the 
source sizes. In particular, the linear priors used bv lBate et al.l (j2008l ) and iFlovd et al.l (j2009l ) will 
favor steep temperature profiles, while a logarithmic prior would favor shallow temperature profiles. 
The same probl em occurs for s parsely sam pled light curves, as illustrated by the studies of X-ray 
microlensing by iMorgan et al.l (l2008bl ) and iDai et al.l (120101 ) , where the lower li mits t o the X-ray 
size a re found to be prior-dependent. From this point of view Bate et al. (j2008l ) and Floyd et al 



(j2009l ) are really upper bounds on /3, and analyses of light curves will be required to determine a 
lower bound. 

Arguments for a flatter emissivity profile in accretion disks have existed for a long time, largely 
based on the mismatch betw een the predi cted and observed spectra of quasars (see the reviews by 
Koratkar Sz Blaes (j 19991 ) and Blaes ( 2004 )). The emission profile can be flattened either by raising 
the temperature of the outer disk (e.g. irradiation of the outer disk by the inner) or by changing the 
balance between radiation and advection in the inner disk (e.g. thick or slim disks). Microlensing 
provides the first probe able to actually measure the physical scales associated with the emission 
regions, and the results are clearly beginning to test the simplest theories. As the numbers of 
systems, the uncertainties in the size measurements and the wavelength range spanned by the 
measurements increases, microlensing can quantitatively address all these issues. 
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Facilities: CTI0:2MASS (ANDICAM), Hiltner (RETROCAM), WIYN (WTTM), HST (NIC- 
MOS, ACS). 
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Table 1. Measured and Derived Quantities 
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103 


2.37 


15.91°:^ 


0.211 


Balmer, Fe II UV 
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Balmer, Fe ii UV 
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Balmer, Fe II Optical 
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0.04 
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Balmer, Fe II UV 


21.97 ±0.19 
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48 
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Note. — Rs from microlensing is the accretion disk size at Ar-est, the rest- frame wavelength corresponding to the center of the monitoring filter 
used for that quasar's light curve. Use half-light radii (J?i/2 = 2.44i?s) to compare these size measurements to other disk models. Significant 
sources of unmicrolensed flux from the QSO Broad Line Region (BLR) that fall into or overlap with the observing pass band are indicated: Balmer 
Continuum (A < 3650A), Fe II UV Continuum (A < 3100A), Fe II Optical Continuum (4240A < A < 5400A) or Mg II (A2798A). Icorr is the corrected 
(unmagnifled) /-band magnitude. Typical /-band measurement errors are < 0.1 mag, but the larger errors on Icorr come from uncertainties in 
the lens magnification. Rs calculated using corrected /-band magnitude and thin disk theory is also unsealed; it is the disk size at the rest-frame 
wavelength corresponding to the center of the H ST I -hand filter (F814W). Both disk sizes assume an average inclination angle i = 60°. 

"The C IV emission line width from lYee fc De Roberti3 l ll99j) depends strongly on the fit to several blended C IV absorption features, so we report 
Mgfi at lower precision. 
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Fig. 1. — Inclination-corrected accretion disk size -R2500 versus black hole mass Mbh- The solid 
line through the data shows our best power-law fit to the data and the dot-dashed line shows the 
prediction from thin disk theory [L/Le = 1 and rj = 0.1). The shaded band surrounding the best 
fit shows the expected variance due to inclination. Disk sizes are corrected to a rest wavelength 
of Xrest = 25OOA and the black hole masses were estimated using emission line widths. The filled 
points without error bars are -R2500 estimates based on the observed, magnification-corrected /-band 
fiuxes. They have typical uncertainties of 0.1-0.2 dex. Solid lines at the bottom of the plot show 
the innermost stable circular orbit for a maximally rotating Kerr black hole and a Schwarzschild 
black hole, representing a plausible range of radii for the inner edge of an accretion disk. 
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Fig. 2. — Results of the power-law fit to i?2500 SlS a function of black hole mass. The contours show 
the 1 — 3a one-parameter confidence intervals for the slope a and the normalization R2bOo{^BH = 
IO^Mq) for the 2500A accretion disk size corresponding to Mbh = lO^M©. The best-fit value is 
indicated with a black point. The filled points along the dot-dashed line are theoretical thin disk 
sizes for quasars radiating at the Eddington limit and with efficiencies of = L/(Mc^) = 0.01, 0.1 
or 1.0. 
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Fig. 3. — Thin disk flux size estimates versus accretion disk sizes from microlensing. For reference, 
the solid line indicates a one-to-one relationship between thin disk Atix size estimates and the 
microlensing measurements. The dot-dashed line is the best fit to the data. Since the data points 
have large errors relative to their dynamic range, the best-fit slope is consistent with unity and its 
average offset from the solid line is 0.6 dex. 
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Fig. 4. — Ratio between the flux and microlensing size estimates as a function of the temperature 
profile slope (3, where T oc . The fiducial power-law model (solid) is normalized by the observed 
ratio at the /3 = 3/4 slope of the thin disk model. The ratio must be shifted into the region delineated 
by the horizontal dashed lines for the two sizes to agree given the uncertainties. Also shown are 
the effects of adding an inner disk edge at Rm = 0-lRx (short dashed curve), making 30% of the 
emission contamination from larger physical scales such as the broad line region (dot-dashed curve) , 
and scattering 30% of the disk fl ux on larger physical sc ales (lo ng dashed curve) . Point s with error 



bars show li mits on from P08 (IPoindexter et al.l 120081'). FOQ (jFlovd et al.l 120091 '). BOB (IBate et al 



20081 ^. AOS (|Anguita et ahlboosl ) and EOS (|Eigenbrod et ahlboosl ) 



